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Abstract 

We consider the effects of neutralino-stau (x — f) coannihilations on the cosmological relic 
density of the lightest supersymmetric particle (LSP) \ in the minimal supersymmetric ex- 
tension of the Standard Model (MSSM), particularly in the constrained MSSM in which 
universal supergravity inputs at the GUT scale are assumed. For much of the parameter 
space in these models, \ is approximately a U(l) gaugino B, and constraints on the cosmo- 
logical relic density Q^h 2 yield an upper bound on mg. We show that in regions of parameter 
space for which the cosmological bound is nearly saturated, coannihilations of the B with 
the f, the next lightest sparticle, are important and may reduce significantly the B relic 
density. Including also B coannihilations with the e and ft, we find that the upper limit on 
m~ is increased from about 200 GeV to about 600 GeV in the constrained MSSM, with a 
similar new upper limit expected in the MSSM. 



Supersymmetry is one of the most appealing options for possible physics beyond the 
Standard Model, motivated theoretically by the help it offers in stabilizing the gauge hier- 
archy, its successful prediction of sm 2 9w in the context of GUTs, and its consistency with 
the range of Higgs boson masses favoured by the precision electroweak data. Accordingly, 
supersymmetry has been the focus of intense phenomenological studies, particularly in the 
framework of the minimal supersymmetric extension of the Standard Model (MSSM) [|TJ. 
Many of these studies assume that supersymmetry is broken by unspecified dynamics in 
some hidden sector of the theory, which is communicated to the observable MSSM particles 
by gravitational interactions. One may further assume that the supersymmetry-breaking 
mass parameters mo, m\/2,A and B are universal at the supergravity input scale, providing 
the constrained MSSM (CMSSM) framework that is privileged in this paper. 

The lightest supersymmetric particle (LSP) is stable in the MSSM, unless additional 
R- violating interactions are postulated 0. The LSP is generally thought to be the lightest 
neutralino x @> U an d is a favoured candidate for the Cold Dark Matter favoured by astro- 
physicists and theorists of structure formation. The relic LSP density can in principle be 
calculated reliably as a function of the parameters of the MSSM [[J. It is remarkable that 
there is a large generic domain of the parameter space in the MSSM and in the CMSSM, 
consistent with all experimental constraints ||, in which the x has a relic mass density 
f2~/i 2 ~ 0.1 as favoured by astrophysical and cosmological arguments. 

The phenomenological arguments for supersymmetry based on the gauge hierarchy, sin 2 8w 
and the Higgs mass nih all suggest that supersymmetric particles should weigh ~ 1 TeV or 
less, but do not provide very precise upper limits on their masses. For example, the amount 
of fine tuning required to maintain the gauge hierarchy increases as the MSSM mass param- 
eters are increased, but there is no objective criterion how much fine tuning is tolerable [0. 
Moreover, sin 2 ^ and rrih are only logarithmically sensitive to the sparticle masses. On 
the other hand, the LSP relic mass density is very sensitive to m~, since the annihilation 
cross section tends to decrease as m~ increases, increasing also its relic number number den- 
sity. The relic density is also very sensitive over much of the parameter space to the scalar 
mass parameters, as these, too, control the annihilation cross section. The constraints on 
the general MSSM parameter space have been explored in some detail, particularly in the 
CMSSM framework. The possibility that the LSP might be mainly a photino 7 has been 
excluded by lower limits on sparticle masses from LEP and elsewhere ||. When studying 
the possibility that the LSP might be largely a Higgsino H, coannihilations || between the 
Higgsino-like LSP and the next-to-lightest supersymmetric particle (NLSP) have to be taken 



into account |9|, |I0 . This Higgsino LSP possibility has also been tightly constrained by LEP 
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and may be explored completely by upcoming runs ||. This leaves us with the likelihood of 
a Bino- (B-) like LSP, and cosmology imposes an important upper limit on its mass, which 



has been given as m B < 300 GeV in the MSSM fTJ and < 200 GeV the CMSSM framework 

u . 

The purpose of this paper is re-evaluate this upper limit, including for the first time 
detailed calculations of coannihilations between the B LSP and the lighter supersymmetric 
partner of the r, the right-handed stau tr, which is the NLSP in much of the B LSP region. 
We find that B — tr coannihilation is particularly important when the relic mass density 
is close to the cosmological upper limit, which we take to be Q~h 2 < 0.3, resulting in a 
considerable relaxation of the previous upper bound on m^. Including also B coannihilation 
with the and ftR, we now find that may be as large as 600 GeV in the CMSSM or 
MSSM. 

As already commented, the LSP is a B in much of parameter space that leads to an 
interesting relic density, both in the generic MSSM and in the CMSSM |J. Indeed, this is a 
prediction of the CMSSM. Unless the B mass happens to lie near m^/2 or m^/2, in which 
case there are large contributions to the annihilation through direct s-channel resonance 
exchange, the dominant contribution to the BB annihilation cross section comes from crossed 
t-channel sfermion exchange. The resonant case is anyway not relevant for the upper bounds 
on rrig to be discussed here. In the absence of such a resonance, the thermally-averaged 
cross section for BB — > ff takes the generic form 
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where Y L ^ is the hypercharge of fh(R), ^/ = m ^ + m \ ~ m % an d we have shown only the 
leading P-wave contribution proportional to x = T/m-g. 

The upper limit on due to the cosmological relic density comes about as follows [|TT|. 
The assumption that the B is the LSP requires, in particular, that < mj. In order to 
minimize the relic density, we must maximize the cross section, which is done by setting 
mj = mg. The cross section is then approximately inversely proportional to m~. The 
cosmological upper limit on fl^h 2 translates into a lower limit on (av) which then, in turn, 
yields an upper limit to mg. In the MSSM, this limit is mg < 300 GeV, when all sfermion 
masses are taken to be equal at the weak scale. 

In the CMSSM, the argument is somewhat similar, although and the sfermion masses 
are no longer entirely independent, because it is assumed in the CMSSM that there is a 
common scalar mass mo at the GUT scale. For a given value of the common gaugino mass 
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my 2 at the GUT scale, the relic B density falls with mo, since m 2 j = Wq + C/m^ 2 + 0(m 
where C/ is a positive numerical coefficient that is calculable via the renormalization-group 
evolution of the sfermion masses. Therefore, the cosmo logical upper limit onQ^h 2 translates 
at fixed my 2 into an upper limit on m . Typically, this upper limit is not larger than 
m < 150 GeV, unless one is sitting on a direct-channel pole, i.e., when mg ~ m^/2 or 
rrih/2, in which case s-channel annihilation is dominant and there is no upper limit to mo- 
However, as already mentioned, this is not our case, as we are interested in an upper bound 
on mg. We recall that scales with my 2 , and it transpires for my 2 ?Z 400 GeV that 
nig exceeds mass of the lightest sfermion, which is typically the tr, for mo small enough to 



satisfy the cosmological bound [12). Thus, the LSP is no longer a neutralino for such large 
values of my 2 , and hence an upper bound < 200 GeV [12| can be established. [] 

When mg attains this upper bound, the B is degenerate in mass with the tr, and 
quite close in mass to the <3r and /2r. It is well known || that, in such circumstances, the 
neutralinos can be maintained in equilibrium by scatterings with a slightly heavier particle, 
and the number density of neutralinos can be significantly reduced by such coannihilations. 
The case of heavy Higgsinos is a well studied example [[H]]. Analogously to that case, the 
B relic density can be reduced through coannihilation with slightly heavier r R 's or other 
sleptons, as we now discuss in detail. 

To derive a thermally-averaged cross section, we use the technique of ||. Thus, we 
expand (<Jv re \) in a Taylor expansion in powers of x = T/m^: 

(av vc i) =a + bx + 0(x 2 ) . (2) 

Repeating the analysis f| for initial particles with different masses mi and m 2 yields 

(^rei) = — ^1 - 2mim2 ) w(s)\ s ^ mi+m2 y +3{mi+m2)T + 0(T ). (3) 

where 



w[s) 



■ \J rfLIPS|M| 2 (4) 

Here dLIPS is the Lorentz-Invariant phase-space element, p(s) is the magnitude of the three 
momentum of one of the initial particles in the center-of-mass frame, as a function of the 



lr This upper bound can be strengthened by requiring that the global minimum of the effective potential 
of the MSSM conserve electric charge and color O] . 
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total center-of-mass energy- squared s, 9qm is the center-of-mass scattering angle, and \M.\ 2 
is the absolute square of the reduced matrix element for the annihilation, summed over final 
spins and averaged over initial spins wherever appropriate. The a and b coefficients in (0) 
may be read off the right-hand side of (^), after expanding in powers of x. f\ 

If the masses of the next-to-lightest sparticles (NLSPs) are close to the LSP mass: 
AM = 0(xf)M, where Xf ~ (1/20 — 1/25) is the value of T/m^ at the time of neutralino 
decoupling, the number densities of the NLSPs have only slight Boltzmann suppressions 
with respect to the LSP number density when the LSP freezes out of chemical equilibrium 
with the thermal bath. In such a case, coannihilations of NLSPs with the LSP, along with 
NLSP-NLSP annihilations, may play an important role in keeping the LSPs in chemical 
equilibrium with the bath ||. These processes can be particularly important when the LSP 
annihilation rate itself is suppressed, as is the case for neutralinos. Gaugino-like neutralinos 
typically annihilate predominantly into fermion pairs, and such processes exhibit P-wave 
suppressions @, so that a <C b in (0). This effect can be seen from ([!]) where the a-term is 
suppressed relative to b by m^/m~, reducing the neutralino annihilation cross section by a 
factor of 0(xf). We also emphasize that 2-2 scatterings with particles in the thermal bath 
keep the NLSPs, in this case the f R , e R and /i R , in chemical equilibrium with each other 
and with x, down to temperatures well below the temperature at which the comoving LSP 
number density freezes out. 

We consider the total density n = J2i n i, where the index i runs over r R , r*, e R , e*, J1 R and 
/I* as well as x, and write the rate equation for n\ 

= ~3Hn - (cr cS v rc i)(n 2 - n 2 eq ), (6) 
where H is the Hubble parameter, and 

^ = Ew- ( 7 ) 

Here = n^ q /n eq where nf 1 is the equilibrium density of particle species i, and is the 
total cross section for particle i to annihilate with particle j. Since the sleptons decay into 
neutralinos after freeze-out, the number density of neutralinos becomes n. Many of the o"^- 
are related, and if we take the e R and J1 R to be degenerate in mass and ignore the electron 
and muon mass, we can write 

CTcff = X X r x r X + 4 a xr r x r r + 8 (J X e r x r e + 2 (°Vr + 0"tt* ) r r?Y + 8 (<7 re + <7 Te * )r r T e + 

4 (cr ee + a ee *)r e r e + 4 (a efl + <r etl * )r e r e (8) 
2 A similar calculation is necessary when the LSP is assumed to be a sneutrino [O . 
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In Table [I] we list the sets of initial and final states for which we compute the annihilation 
cross sections. The cross sections for other reactions, such as £ R £* — > hH, hA, hZ, H + H , 
W + H~, . . . are either suppressed or kinematically unavailable in the regions of CMSSM 
parameter space relevant to our analysis. In practice, we find that the dominant contributions 
to cr eff come from annihilations of £^£^ to gauge bosons, i^i^ to lepton pairs, and £^x t° 
t + gauge boson. 

Table 1: Initial and Final States for Coannihilation: {i,j = r,e,fi} 



Initial State 


Final States 


pi Ji* 
R l R 


77, ZZ, iZ, W + W~, hh, tt 






R *-R 




^R ^R 1 * 7^ 3 




Ax 


t-i.tZ.th 



To get a simple estimate for the size of the effect of including the next-to-lightest states, 
we first assume degenerate LSP and NLSPs, and consider a model in which the NLSP-NLSP 
and NLSP-LSP annihilations are all unsuppressed. Thus, we take 

{dij ps aij ^ feu, an = 0;i,j > 1}, (9) 

where the subscripts i,j refer to the NLSPs and 1 to the LSP. Denoting with superscripts 
quantities that are computed ignoring the NLSP states, we estimate the following ratio of 
relic densities without and with coannihilation: 

^fr^O'l'- (10) 

where a$/x f ps 1 + x^\n(g tot / g^) w 1.2, g tot = Y,i9h and a eff /feii w 1 - gl/sht = 15/16 
for the case of three degenerate slepton NLSPs. Thus, in this crude approximation we find 
a factor ~ 35 reduction in the relic density. Ignoring the (heavier) left-handed sleptons, 
we may reduce ([[]) to (av) ~ 3gfx/(16iim^), yielding fi^h 2 ~ 8 x 1CT 6 m|. Thus, in this 
simple approximation, Q^h 2 = (Q~h 2 )°/R < 0.3 gives an upper bound on the B mass of 
m~ < 1.2 TeV. 

We have gone beyond the above crude approximations to make a detailed numerical 
analysis of coannihilation effects on the neutralino relic density, including light sleptons, 



5 



some of whose results are displayed in Fig. [I]. The light shaded region corresponds to 
0.1 < ^^h 2 < 0.3, and the dark shaded region to < m%. We have chosen the repre- 
sentative points tan/3 = 3 and 10, and present results for both \x < and fi > 0. In practice, 
the relationship (Q) is not exact, not all of the are unsuppressed, there are contributions 
from the b^, and the e R and fl R can be slightly heavier than the f R . These corrections have 
the net effect of reducing a e g by a factor of ~ (3 — 4). Numerically, we find that R ~ 10 along 
the line where = m~ , at the top of the dark shaded region in Fig. |. As m increases, 
the sleptons become heavier relative to the neutralino, and their number density rapidly 
falls, reducing their contribution to a e g. The relic density rises rapidly in this region, leaving 
an allowed band in mo which is about 30-50 GeV wide for < 800 GeV. In Fig. ||] we 
extend the coverage of Fig. [l|a,c over a larger scale, to show the cross-over point between the 
regions with ^h 2 < 0.3 and m~ R < m%, where there would be an unacceptable abundance 
of charged dark matter |Q. The two constraints together require my 2 ^ 1450, corresponding 
to an upper bound on the neutralino mass of m^ < 600 GeV. The results for fi > and 
fi < are very similar, so we do not display the latter. The width of the allowed region is 
insensitive to Ao, though the position of the line m^ = can vary somewhat. As already 
commented, the requirement that the electroweak vacuum conserve electric charge and color 
constrains significantly the CMSSM parameter space [fT3|| . We find that the large-mi/2 tail 
of the region newly allowed by coannihilation obeys this requirement, for tan/5 > 3 for some 
values of A . 

We expect the corresponding bound in the MSSM to be very similar. In the general case, 
one must take all the squarks and sleptons degenerate with the neutralino and compute 
the annihilation and coannihilation cross-sections for all possible combination of sfermions. 
However, if the rates are the same as for the sleptons, the effect is about a 15% decrease in 
(Q~h 2 )°/R, leading to a similar bound on m^ as in the CMSSM. 

The potential significance of coannihilation effects had been emphasized previously || 
H |i~0 |, particularly in the Higgsino LSP region. We have shown that this can also be an 
important effect in the B region, where the NLSPs include the tr and other right-handed 
sleptons. In particular, we find in the context of the CMSSM that the upper bound on the 
LSP mass quoted previously [0 should be increased by a factor of about two, to < 
600 GeV. 

This observation has many potential ramifications, in particular for searches for neu- 
tralino dark matter. Generically, in regions of parameter space where coannihilation is 
important, a relic density of astrophysical interest is now obtained for a smaller annihilation 
cross section. This is likely to reduce the typical rates for signatures of annihilations in the 
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galactic halo. The corresponding elastic scattering cross section is also likely to be reduced 
generically, with a consequent suppression of signatures for scattering on nuclei and capture 
followed by annihilation in the Sun or Earth. These points are worthy of further study. 

On the basis of the previously-quoted upper bound on the LSP mass, it has been argued 
that the LHC is guaranteed to detect some supersymmetric particles |{L5|] , since it could reach 
out to mi/2 ~ 1200 GeV when mg ~ rrif R . It is still true that the overwhelming majority of 
the CMSSM parameter space allowed by cosmology can be explored by the LHC, but there 
may be a narrow region of mo extending up to my 2 ~ 1500 GeV that is problematical. This 
point is also worthy of further study. 
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Figure 1: The light-shaded area is the cosmologically preferred region with 0.1 < < 0.3. 

The dashed line shows the location of the cosmologically preferred region if one ignores the 
light sleptons. In the dark shaded region in the bottom right of each panel, the LSP is the 
f R , leading to an unacceptable abundance of charged dark matter. Also shown as a dotted 



line is the contour m~ = 1.1 m^. 
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